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ON  THE  ITERATION  OF  ?a;EK  SERIEO  IN  TVO  VARIABLES 
Richard  Bellman 

§1.  Introduction 

It  was  shown  by  Koenigs,  C  5^  ,  that  if 

(1)  f(x)  -  X 

k-1 

has  a  non-zero  radius  of  convergence  and  |a^|  is  leas  than  1  and  not  equal  to 
zero,  then  for  |x|  sufficiently  small,  there  exists  a  function  of  two  variables, 
f(x,  t)  possessing  the  properties  of  a  generalized  iterate,  which  Is  to  say 

(2)  f(x,  0)  -  X, 

f(x,  n)  •  n-th  iterate  of  x,  for  n  a  positive  integer, 
f(f(x,  s),  t)  -  f(x,  s  ♦  t),  fors,  t>0. 

This  function  has  the  elegant  representation, 

(3)  f(*,  t)  -  ^^|\^^(x)J, 

where 

U)  ^(x)  -  lia  f^"^(x)/^. 

n  -eoo 

The  function  J^(x)  is  Itself  an  interesting  and  important  function,  since  it 

■linearizes*  f(x),  i.e., 

(5)  ^(f(x))  -  • 

The  problem  of  finding  a  ^(x)  satisfying  (5),  given  f(x)  has  a  long  history, 
going  beck  to  Abel.  For  further  discussion  and  references,  we  refer  to 

Hedamard,  [4]  . 
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Let  US  now  turn  our  attention  to  power  series  in  two  variables.  Let 


(6) 


f(x.  k,  />0,  kW>l 

g(x,  y)  -  X  > 


be  convergent,  for  |x|  and  |y|  sufficiently  small  and  let  us  make  the  assumption 
itiit  the  characteristic  roots  of  the  matrix  of  the  coefficients  of  the  linear 
terns, 


are  non-2eru  and  less  than  unity  in  absolute  value. 

The  iterates  of  the  function-pair  (f,  g)  are  defined  as  follows: 


(^)  *  ^n-1^^’ 

gn(x,  y)  -  6^1-1^^’  S),  n  -  1,  2,  •••, 

where  ve  set  y)  -  x,  gQ(x,  y)  -  y. 

It  is  natur.l  to  ask  whether  we  can  find  continuous  iterates,  i.e,,  a  pair 
of  fur.ctions  of  three  variables,  ^f(x,  y,  t),  g(x,  y,  t)^,  for  which 


(9) 

f(x,  y 

»  ■ 

ryx, 

y) 

g('<i  y 

,  - 

S„(x, 

y). 

n  '  0,  1, 

i  •  •  • 

and 

(10) 

f  f(x, 

y»  3)i 

g(x, 

y,  s),  tl 

-  f(x,  y,  3 

♦  t). 

g'  f(x, 

y,  s), 

g(x, 

y »  »  i ) 

•  g(x.  y»  9 

♦  i), 

for  s  an  a  t  >0. 


In  general,  t.hese  functions  will  be  defined  only  for  x  and 


Ixl  sufficientlj  small  and  s,  t  >  0. 

Following  Koenigs'  idea,  which  is  really  an  application  of  the  method  of 
successire  approximations,  the  existence  of  these  functions,  together  with  an 
•xi<licit  representation,  will  ♦‘ollow  readily  i  f  we  establish  the  existence  of 
functions,  ^(x,  y),  4«(x,  y),  satisfying  the  malogue  of  (5),  namely,  the 
rector-matrix  equation 


Before  proceeding  any  f’virther,  lot  us  note  that  by  means  of  some  elementary 
transformations,  wo  may  reduce  the  problem  to  the  case  where  f  and  g  have  the 

sls^xler  forms, 

(12)  f  -  px  ♦  X 

k*/>2 

rg  -  oy  ♦  X  »  (i,  o  /  0, 

k*£>2 

prorldeo  that  A  does  not  have  multiple  characteristic  roots.  If  A  has  multiple 
characteristic  roots  a  representation  similar  to  that  given  above  holds  and  the 

details  are  very  much  the  same. 

Under  the  assiaption  that  A  has  distinct  characteristic  roots,  the  equations 
in  (U)  take  the  simpler  form 

(13)  g)  -  Pl^(x,  y), 

♦(r»  f)  -  o4-(x»  y). 

By  analogy  with  (4),  we  might  expect  that 
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(14)  ^  -  lljn  y)  /  p" 

n  ->oo 

4  -  llin  g^(x,  y)  /o"  . 

n  ^00 

Unfortunately,  very  sinple  examples  show  that  these  limits  may  not  exist.  For 
example,  consider 

(15)  f  -  px 

g  •  ay  *  . 

It  is  easily  shown  by  induction  that 

(16)  f^  -  p"x 

n  /  n-1  n-2  2  2n-.  2 

-  0  y  ♦  (o  ♦  0  p  ♦  •••  ♦  p  )x  . 

Thus , 

(!'!’)  ^  “  lim  /  p”  ■  X 

n  ^00 

2  2 

ip  ■  11m  ■  y  ♦  “  lim  1  ♦  ^  ♦  •••  ♦ 

0  _  0 

n  ->CD  n  -►CD 

This  last  limit  exists  only  for  \pi'o\  <  1. 

Despite  the  fact  that  a  direct  application  of  Koenigs'  method  falls,  proofs 
of  the  existence  of  functions  satisfying  (13)  have  been  given,  provided  that  a 
certain  condition  is  satisfied  by  p  ani  a  whicr.  we  shall  discuss  below,  using  the 
method  of  successive  approximations,  or  the  method  of  majorants,  the  last  being 
particularly  suited  to  the  analytic  case  wy.ich  we  are  treating  here,  cf. 

Leau ,  [  6]  ,  Drevy,  [  3  . 
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If  w«  Miploy  the  method  of  undetermined  coefficients  to  find  functions  ^ 
•nd  ^  satisfying  (13)  for  the  functions  f  and  g  defined  by  (15),  we  find  without 
difficulty, 

(18)  i  -  X 


4-  -  y 


If  "  V 


o  -  p 


m 


provided  that  a  ^  p  , 

If  we  compare  the  solution  of  (13)  with  the  att«if)ted  solution  of  (17),  we 
are  at  once  struck  by  the  fact  that  1/  (o  -  p^)  is  the  generalized  sum  of 

OD  2  n 

X  (P  /o)  /  o  ,  using  an  appropriate  sxaramability  method,  or,  alternatively, 
n  *0 

00 

1/  (1  -  z)  is  the  analytic  continuation  of  the  f\inction  defined  by  ^  z”  for 

n«0 

l»l  <  1. 

This  fact  leads  us  to  conjecture  that  a  suitable  interpretation  of  the  limits 
in  (17)  will  yield  the  required  functions  i  and  4.,  and  we  shall  show  below  that 
this  is  actually  so.  Although  the  method  we  employ  Is  equally  applicable  to  the 
general  case  of  power  series  in  n  variables,  for  the  sake  of  simplicity  we  restrict 
ourselves  to  the  two-variable  case. 

We  feel  that  this  method  of  establishing  the  existence  of  the  linearizing 
functions  4  vid  is  of  interest  because  it  furnishes  another  example  of  the 
Mta-aathematical  axiom  that  formal  methods  correctly  interpreted  yield  correct 
rssults.  An  example  was  given  in  a  pz^vlous  paper,  [l],  where  it  was  shown  that 
the  Liouville-*Neumann  solution 

t  •  g  *  ^J'  K(x,  y)g(y)dy  ♦  ••• 


of  the  Fredholn  integral  equation 


(-^0) 


f  ■  g(x)  ♦  ^  f  K(x,  y)f(y)dy 
J  0 


is  suninable  by  an  appropriate  aumability  method  to  the  Fredholm  solution 
whenever  X.  is  not  equal  to  a  characteristic  value. 

In  many  cases,  these  formal  methods,  combined  with  sunaability  methods, 
are  useful  for  computational  purposes,  cf.  Buchner,  [2], 


^2,  Some  Intermediary  Functions 
Let  us  set 


(1)  f  -  px(l  ♦  u)  ♦  r(y) 

g  -  oy(l  ♦  v)  ♦  8(x), 

where  u  and  v  are  power  series  in  x  and  y  lacking  constant  terms  and  r  and  ■ 

are  power  series  in  y  and  x,  respectively ,  lacking  constant  and  first  degree 

terms.  We  shall  employ  the  following  notation  for  the  iterates: 

fn.l  -  «n)'  Vl  '  «n)> 

•n.l  ■  '•n.l  • 

Vl  ■  «n>-  Vl  ■  «n)>  "  '  1. 

where  f^  ■  ■  g.  Using  this  notation,  we  have,  for  example, 

2 

f^  -  p  x(l  ♦  u^){l  ♦  u^)  ♦  pr^(l  ♦  u^)  ♦  r^, 

g2  -  o^yd  ♦  v^)(l  ♦  v^)  ♦  08^(1  ♦  T^)  ♦  s^, 


(3) 


and  it  is  not  difficult  to  establish  inductively  that 


ik) 


-  p”x  IT  (1  ♦  T  ♦  \) 

k"l  k-2 


♦  • . .  ♦  r 


n 


or 


(5) 


n 


T  (1  •  %) 

k-l  “ 


n 

/-I 


T  (1  ♦  \) 

k  -  i 


and,  similarly. 


(6) 


h. 


00 

0  k-l 


(1  ♦ »,) 


n 

y  -r 


o'^T  (1  •  vj 

k-l 


If  whs  infinite  products  converge,  and  also  the  infinite  series,  we  may 
pass  to  ths  limit  as  n— ^oo  and  obtain  i  and  a..  It  is  easy  to  demonstrate, 
aftar  the  fashion  of  Koenigs,  that  for  |xj  and  |y|  sufficie:*!/  small,  the 
infinite  products  converge,  on  the  assumption  that  '^|,  |o|  are  both  less  than 
one.  He  have  already  assumed  thit  p  ^  o,  and  without  loss  of  renerality  we  may 
take  |p|  >  |o|.  Actually,  in  place  of  p  /  o,  one  need  jnly  assume  th=it  the 

■atrix  A  is  diagonalizable. 

2 

|P|  't  |o|  >  IPI  »  infinite  series  which  occur  in  the  limiting  forms 
of  (5)  and  (6)  converge.  If  o  ■  p'^,  n  -  3,  **•,  the  functions  4  an!  ij  need 

not  axist,  as  we  see  by  taking 

(7)  f  -  px 

00 

g  -  oy  ♦  2  c  x*", 
k  -  2  ^ 


for  which  4  and  are  given  by 


(8) 


4  - 


CD  C,  X 

t.  -y  •  Z 


k  -  2  p  -  0 


2  2  3  . 

To  shorten  our  proof,  we  shall  assume  a  4  9  »  \9\  >  |ol  ^  iPl  »  '^hich  is  the 

first  interesting  C4se, 

Let  us  define 


CD 


(9) 


F(z»  y)  -  2! 


n  n 


n  ■! 


n 


T  (1  ♦ 

k  -1 


n 


00  8  S 

G(zi  X,  y)  •  X  -T - - - 

TT  (1  ^  V,) 

k-l 


The  convergence  of  these  functions  for  |x!,  |y|  small  and  |z|  <  1/  |pj‘ 

follows  from  inequalities  for  r  ana  s  we  derive  below.  Vie  also  require  the 

n  n 

functions 


00 


(10) 


itj  y)  •  T  (1  •  “k)[*  ' 
k-l 


CD 


4.(2,  X,  y)  -  TT  (1  ♦  ^ 

k-l 


Lenuna  1 .  For  | z  1  <  1  /  •  p l  , 


00 


4{-,  f|  g)  ■  J  4iz,  X,  y)  ♦  x(p  -  ■j)  TT  (1  ♦  Uj^) 

k  -1 


\  1  ® 

4'(z.  f»  g)  -  7  y)  ^  -  7)  TT  (1  ♦  yJ  . 


k-l 


(11) 


225 


-9- 


Proof  of  Lenut:  We  have 


(12) 


00 


F(z,  g) 


n 


y  — ! 

ir  (1  •  V 


(1  •  u^) 


F(z>  y)  - 


k-2 


whence 

(13) 


00 


C)  •  TT  (1  ♦  u,  )  [f  ♦  F(z,  f,  g)  , 
k-2 


f  (1  ♦  u,  ) 

k-1 


^  I 

j  ,  T  ♦  -  F(z,  X,  y) 

V  1  ♦  )  2  '  '  *  ' 


f  (i  ^  U,) 

k-1  ^ 


1  i 

px  -  -I 


/ 


f  (1  *  U^) 

k  -1 


-  X 


^  ♦  X{p  -  i)  ff  (1  ♦  U^), 


k-1 


and  similarly  for  <f. 

From  (11),  we  see  thit  if  4-  possesses  in  'in'iiytic  continuation  which 
includes  the  point  z  •  1/a,  then  we  will  have 


(14) 


^(“.  «.  y)  •  p^(“.  'i  y'h 

P  P 

4’(^»  X,  y)  -  o«p(^,  X,  y), 


and  consequently  the  desired  linearizing  functions. 

In  the  next  section  we  prove  that  this  ir.ilytic  continuation  exists. 


§3#  Proof  of  Analytic  Continuation 
Wa  wish  to  demor.strate: 
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Loaaa  2.  The  functions  1*^(2,  x,  y),  y)»  considered  as  functions  of  s 

for  I X '  and  1  y |  sufficiently  small  are  meronorphic  functions  whose  only 

possible  Si.  r^ularities  are  at  the  points  z-o  'p  ,  m.n>0  for 

.  -(n>2)  HB  ^  ri  r  . 

and  2  -  p  'a  ,  m,  n  >  0  for  4<. 

The  proof  depends  upon  a  continued  application  of  the  Hadamard  multiplication 
theorem  which  we  state  as 

Lemma  3.  (Hadamard) .  IJJ  ^  1 2  I  1 ,  the  3in;:ularitie3  of  the  functions  defined 
00  00 

^  f(z)  -  2  a  gizj  -  2  ^^2",  for  Izi  <  a  ^  1,  are  at  a . |  3  , »  respectively, 

n-0  "  n-0  ^ 

1  -  1,  2,  •••,  k,  J  -  1,  2,  •••,  /,  then  the  singularities  of  t he  function 

00 

define  d  ^  hiz)  -  2!  a  for  |z|  a  <  1,  are  to  be  found  in  the  region  |z|  1 

at  the  points  if  any  singularities  are  present, 

^  %i 

We  require  first  the  following  crude  inequalities: 


Knl 

<  pjl 

X| 

«rl 

<  ojl 

for  1  y  1  < 

Ix| 

i  s. 

'fpl 

<  Pi' 

y  1 

V  0") 

y  f 

for  *  X  1  ^ 

'y  1 

<  s. 

where  P^  ■  1  p  1  ♦  ^  1 

" 

1  0  1  ♦  £ 

and  can 

be 

made  arbitrarily  small  by  taking 

S  sufficiently  small.  These  inequalities  follow  from 

(2)  |f|£  p!|x|(l*iui)^aj^yi^ 

<  Ip  I  'X  I  '^1  ♦  |u  !  ♦  a^lyKi  y'  /  |x|)^]  <  p^iX  I 
if  y'  ^  lx  ^  §.  Similarly,  we  obtain 
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(3)  |g|<  io||y|(l  ♦  'vl)  oilxlj  -  lv')  ♦bj^ixl 

<  0^1 X  !  . 

From  these  results,  (1)  follows  ty  iteration. 

As  aientioned  above,  we  shall  consider  only  the  '' ise  where  |^|^  •  o  ^  I  1  ^  • 

The  series  for  F  in  (9)  of  §<2,  is  majorized,  Lecause  of  tne  convergence  of  the 

infinite  product  by 

(O  Fj  -  f  irji"  , 

n  -j. 

which  since  |r^|  <  <  c^o^^lxl,  ( assujuin,.’,  |Xj  >  [y\,,  is  converf;ent  for 

|s|  <1/0^^  and  hence  for  z  -  1/^  if  c  is  small  enouj^h.  The  investi (Ration  of 
0{»)  Is  more  difficult.  Since  3(x)  ■  e  ♦  e_x^  ♦  we  have  s  ■  e,f^  ♦  ^(If  1^) 

•  *2^  *  Hence,  convert:es  for  z\  <1/  anl  thus  for  2  "  . 

if  6  Is  small  enough.  To  ascertain  the  anal/tlc  tehavior  of  G  at  2  ■  1 ''0,  it 
is  sufficient  then  to  consider  the  series 

(5)  0^  -  f  ^"/1r  (1  .  u,y  . 

n-1  K  -1 

In  rlew  of  the  Hadamard  multipl icatlcn  tneorera,  it  is  "of fi'-ient  to  examir.*-  the 

two  series 

00  n  OD  n 

(6)  a,  ■  Z  TT  (1  •  --  -  Z  t:'’'  T  . 

n^lk*l  n-i.  K  ■! 


Dsing  (5)  of  5^,  we  obtain  for  G^  the  simpler  series 
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It  follow^i  that  Qj  has  a  sinple  pole  at  z  •  1/p  and  no  other  singiilarity  within 


<  .  Turning  now  to 


OD 


(8) 


f  (1  ^  V 

k  -1 


00 

^1  ® 

z 

«7Z 

(1 

•  V 

n»l 

!  k-r*l 

00 

00 

^  1’  - 

Z 

n^l 

“n.l 

I 


where  o 


n 


Hence 


(9) 


"(.I. '(5 

OO  y  OD  V  00 

T  (i-uj^)(i-Z\)  IT 

-  k»l  '  1  '  k«l 


00 


1  -  z 


1-1 


SI,  V"  •  - 

n-1 


where  has  no  slng\ilarity  within  |zl  <  l/pj^.  Hence,  we  see  that  ^2  * 

singularity  at  z  -.I  aiid  no  other  within  l/pj^.  From  the  representation 

00  /  f 

(10)  0^  -  E  i  T— -  T  (1  .  V*  • 

"•i\  T  (1  • 

k-1  / 

2 

and  Leaca  3>  i*  follows  that  has  a  possible  singularity  at  z  ■  1/p  and  no 
3  2 

other  within  l/p^.  That  this  singularity  at  z  -  l/p  follows  from  the  fact  that 

2 

the  raaius  of  convergence  is  z  -  l/p  . 

The  merooorphic  behavior  over  the  whole  plane  Is  obtained  similarly,  using 
induction. 
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